Inverse Design of Low-Boom Supersonic Concepts Using 
Reversed Equivalent-Area Targets* 


WuLi 1 

NASA Langley Research Center, Hampton, Virginia 23681, USA 
Sriram Rallabhandi 2 

National Institute of Aerospace, Hampton, Virginia 23666, USA 


A promising path for developing a low-boom configuration is a multifidelity approach 
that (1) starts from a low-fidelity low-boom design, (2) refines the low-fidelity design with 
computational fluid dynamics (CFD) equivalent-area (A e ) analysis, and (3) improves the 
design with sonic-boom analysis by using CFD off-body pressure distributions. The focus of 
this paper is on the third step of this approach, in which the design is improved with sonic - 
boom analysis through the use of CFD calculations. A new inverse design process for off- 
body pressure tailoring is formulated and demonstrated with a low -boom supersonic 
configuration that was developed by using the mixed-fidelity design method with CFD A e 
analysis. The new inverse design process uses the reverse propagation of the pressure 
distribution (dp/p) from a mid-field location to a near-field location, converts the near-field 
dp/p into an equivalent-area distribution, generates a low-boom target for the reversed 
equivalent area (A e?r ) of the configuration, and modifies the configuration to minimize the 
differences between the configuration’s A e?r and the low-boom target. The new inverse 
design process is used to modify a supersonic demonstrator concept for a cruise Mach 
number of 1.6 and a cruise weight of 30,000 lb. The modified configuration has a fully 
shaped ground signature that has a perceived loudness (PLdB) value of 78.5, while the 
original configuration has a partially shaped aft signature with a PLdB of 82.3. 
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Nomenclature 

amplitude of oscillation 

classical equivalent area defined by Mach angle cut method 

equivalent area defined by using reverse propagation of off -body pressure 

body length of the configuration or the aircraft concept 

dimensionless dispersion parameter 

ambient speed of sound 

Mach number 

(calculated pressure - free-stream pres sure)/(free- stream pressure) 
nondimensional pressure or waveform 
ray tube areas 

coordinates of a point on the plane 
relaxation factor, 0.5 

positive constant for regularization of diffusion equation 
dimensionless thermoviscous parameter 
ambient atmospheric density 
nondimensional distance during propagation 
nondimensional time coordinate of the waveform 
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a constant gap between two nondimensional time coordinates 
dimensionless relaxation time parameter 


I. Introduction 

D EVELOPING a configuration with a shaped ground signature that is propagated from computational fluid 
dynamics (CFD) off-body pressure distributions is still a difficult task. A promising path for developing a low- 
boom configuration is a multifidelity approach that (1) starts from a low-fidelity, low-boom design, (2) refines the 
low-fidelity design with CFD equivalent-area (A e ) analysis, and (3) improves the design with sonic -boom analysis 
by using CFD off-body pressure distributions. Steps (1) and (2), along with the analysis portion of step (3), have 
been implemented in ModelCenter 1 as a mixed-fidelity design process for low-boom supersonic concepts. 2 Case 
studies using this mixed-fidelity process have been documented in Refs. 2 and 3. The main focus of this paper is on 
the development of an inverse design process for step (3) of the multifidelity design approach. A case study 
demonstrates how to tailor the off-body pressure to obtain a low-boom configuration with a fully shaped ground 
signature. 

The theory behind the new off-body pressure tailoring is based on a reversed equivalent area (A e?r ) that 
represents an accurate far-field approximation of the aircraft as a body of revolution. The A e?r is calculated from a 
mid-field pressure distribution in the following manner: the mid-field pressure distribution is reversely propagated to 
a near- field location by using the augmented Burgers equation and then converted to the equivalent-area distribution 
by using classical far-field theory. By using A e?r , one can convert off-body pressure target matching to A e?r target 
matching. The new inverse design process starts with calculation of A e?r for a given configuration, generates a low- 
boom target for the configuration’s A e , r , and modifies the shape of the configuration to minimize the differences 
between the configuration’s A e?r and the low-boom target. 

A design case is used to demonstrate that the A e?r target matching problem can be solved with volume A e shaping 
and shape optimization for lift tailoring with three twist angles of the horizontal tail. In this paper, the mid-field 
location is 3BL (or 381 ft) below the aircraft, and the near-field location is 50 ft below the aircraft. The baseline 
configuration was originally designed by using the mixed-fidelity design process 2 for a cruise weight of 30,000 lb, a 
cruise Mach number of 1.6, and a cruise altitude of 45,000 ft. Sonic-boom analysis with CFD off-body pressure 
distributions was used to verify that this mixed-fidelity design had a ground signature with a perceived loudness of 
82.3 PLdB. After the A e?r target matching, the ground signature of the modified configuration is fully shaped and has 
a perceived loudness of 78.5 PLdB. 

In the literature, there were many papers on design methods for developing a supersonic configuration with a 
shaped ground signature. See Ref. 4 for a survey and Refs. 5-10 for some recent advances on the development of 
low-boom design methods. So far, there is no published reference on any demonstration of using a design method to 
generate a low-boom supersonic configuration that has a fully shaped ground signature. 

This paper is organized as follows. Section II introduces the augmented Burgers equation for reverse 
propagation, as well as its solution. In section III, the relationships among A e?r , A e , and the mid-field pressure 
distribution are quantified by examples. The design case is documented in section IV, and concluding remarks are 
given in section V. 


II. Augmented Burgers Equation for Reverse Propagation 

The augmented Burgers equation is used for propagation of the off-body pressure distributions, as well as the 
equivalent-area or F - function distributions. 11 This section describes the augmented Burgers equation and its solution 
method for reversely propagating the mid-field pressure distributions to a near-field location. The reverse 
propagation method yields a more accurate equivalent-area approximation of an aircraft for low-boom design than 
the traditional A e (see section III). 

A. Augmented Burgers Equation 

The augmented Burgers equation includes effects such as nonlinearity, thermoviscous absorption, and molecular 
relaxation phenomena during the propagation of pressure distributions through the atmosphere. The dimensionless 
Burgers equation is shown in Eq. (1), where the last two terms represent geometrical spreading and propagation 
through a stratified atmosphere. 
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As the initial waveform, for which the nondimen sional version is P(0,r), moves downward toward the ground, 
the nondimensional distance a from the waveform to the aircraft increases. The waveform is nondimensionalized by 
using a reference pressure P re f . The shock formation distance of a plane wave with amplitude P re f is used to 
nondimensionalize the actual distance. The dimensionless time r is a function of the arrival time of the waveform, at 
a fixed altitude, for a stationary observer. The dimensionless thermoviscous parameter r, the dimensionless 
relaxation time parameter 6 V , the dimensionless dispersion parameter C v , the ray tube area S , the atmospheric density 
Po, and the speed of sound c 0 are obtained under the assumption of a standard atmosphere. A ray tracing model, 
which computes the ray tube area as well as the atmosphere-related parameters such as molecular relaxation and 
absorption values, is included in the propagation process. The initial waveform is then propagated toward the ground 
or any mid-field location iteratively for a finite number of shock formation distances: 0 = a 0 < o x < ••• < o m = the 
ground or the mid-field location. 

The augmented Burgers equation can be modified for the reverse propagation of the waveform as 
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In contrast to Eq. (1), o = 0 in Eq. (2) corresponds to the waveform at the initial mid-field location, and as a 
increases, the waveform moves toward the aircraft. 


B. Numerical Solution of Augmented Burgers Equation for Reverse Propagation 

While nonlinear effects during propagation steepen the waveform and impose an upper bound on the allowed 
spatial step size during iterations, reversing the nonlinear propagation smoothens the waveform and poses no 
numerical difficulty or spatial step size restriction. The absorption and the relaxation in the reverse direction increase 
the magnitude of the waveform rather than attenuate it, as would occur when these are applied in the forward 
propagation. The same ray tracing routine for the forward propagation is used for the reverse propagation, except 
that the direction of propagation is reversed. 

As in the case of the forward propagation process, the reverse propagation is solved by using an operator 
splitting technique. If both the time-step size Arand the spatial step size (a k -a k - i) are small enough, then the error 
that is induced by such a splitting methodology is insignificant. The reverse nonlinear equation is 
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Equation (3) is solved with an upwind flux-splitting method that uses the following formulas: 
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In the above equations, the symbol cr k denotes the next altitude location in reverse propagation, Aris the constant 
distance between two dimensionless time values for all cr, = P(<J k -i , /• At), and P^u ~ P(o^, i- A t). Note that Eq. 
(4) computes P^ X i explicitly from the known data. 
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Absorption is modeled as a diffusion equation, while relaxation is modeled as a diffusion equation with a mixed- 
derivative term. The reverse diffusion equation is inherently ill-posed, which makes Eq. (2) difficult to solve 
numerically. Several regularization techniques can be used to develop a reverse propagation algorithm that includes 
absorption and relaxation effects. We studied the pseudo-parabolic and quasi-reversible regularization schemes 12 
that are given in Eqs. (5) and (6), respectively. 


dP _ 
~do ~ 

-A — -e 

A dT 2 £ 

■A^( d -) 
dr 2 \doJ 

(5) 

dP _ 

A d2p F 

■ A 92 f d2p ) 

(6) 

do 

A dr 2 £ 

dT 2 \dr 2 J 


In the above two equations, the original diffusion equation is modified by the last term to stabilize the numerical 
solution of the modified equation at the expense of adding numerical error in the solution. The constant s is a 
positive scalar that makes Eqs. (5) and (6) numerically solvable. A desirable value of £ (i.e., 10” 4 for the reverse 
propagation calculation in this paper) is determined by a trial-and-error method: s must be large enough that the 
numerical solution of the modified equation can be generated, and s must be small enough that the solution of the 
modified equation is an accurate approximation of the original diffusion equation. During the course of the 
evaluation and verification, the pseudo-parabolic scheme was shown to be better than the quasi-reversible scheme in 
terms of stability and accuracy. Therefore, all of the reverse propagation results in this paper are generated with the 
pseudo-parabolic scheme. 

The vector P k \ that is defined by Eq. (4) is a solution of the reversed nonlinear equation and can be updated by 
using the following linear equations to include the absorption and relaxation effects: 
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Equation (7) represents the reversed molecular relaxation by oxygen, Eq. (8) represents the reversed molecular 
relaxation by nitrogen, and Eq. (9) represents the reversed thermo viscous absorption. The parameters # v?ox , Q ox , 6t >ni , 
andQ n i are dimensionless relaxation time parameters and dispersion parameters for oxygen and nitrogen, 
respectively. The tri diagonal matrices Pmr(^v,ox?^v,ox)> ^v,m), Pmr(^v,ni?^v,ni)> At vr , and Ptvr *u*e 

defined in the appendix. The final iterative waveform solution P k for the reverse propagation at the altitude location 
c 7 k is given in Eq. (10). The constant >^ st is determined by using the ray tracing method and accounts for the ray tube 
spreading and stratification in the reverse direction. 

The reverse propagation process starts with the waveform P f = P( 0, i- A t) at the initial mid-field location; then 
the waveform P k is computed at the altitude location <j k via the intermediate waveform solutions P k h P k x , P k \, and 
P t y r , for k = 1,2, .... The reverse propagation process terminates when <j k reaches the desired near-field location. 

III. Equivalent Area Defined by Reverse Propagation 

By using the reverse propagation method that is discussed in the previous section, a pressure distribution at a 
mid-field location can be reversely propagated to a near-field location. The near-field dp/p can be converted back to 
an equivalent-area distribution, which is called the reversed equivalent area A e , r , by using the following equation: 

i 

Ae.rW = 4 ' 12 Jo (dp / pW - *o) ■ V* - 1 dt (11) 

In Eq. (1 1), M is the cruise Mach number, R is the distance from the near- field location of the dp/p to the aircraft, x 0 
is the starting location of the nonzero segment of the near- field dp/p, and y — 1.4. Equation (11) is a composition of 
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the relationship between the F-function and dp/p (see Eq. (1) in Ref. 13) and the conversion from F - function to 
equivalent area (see Eq. (3) in Ref. 14). In this paper, A e?r is always defined by Eq. (11) and the reverse propagation 
of the off-body pressure at 3BL below the aircraft to the near-field location R = 50 ft. 

In traditional far-field theory, A e accounts only for the volume (monopole) and lift (dipole) effects, with the 
assumption that the higher order effects are insignificant for sonic -boom analysis. However, the far-field theory 
assumption is invalid for complex, non-slender configurations. The following example quantifies the inaccuracy of 
this assumption for a wing-body configuration. 

Figure 1 shows that the body of revolution that is defined by A e generates a dp/p at 3BL that is significantly 
different from the dp/p of the three-dimensional configuration for 80 < v < 100. The two signatures for the CFD A e 
and dp/p, respectively, are also different, as shown in Fig. 2. The dp/p at 3BL below the configuration is computed 
by using Cart3D 15 with a rotated and stretched volume grid of 20 million cells. 16 

In contrast, the reversed equivalent area is a highly accurate representation of the three-dimensional 
configuration for sonic -boom analysis. Figures 3 and 4 show that A e?r has an almost identical dp/p at 3BF and the 
same ground signature as the wing-body configuration when the CFD off-body dp/p distribution is used for sonic- 
boom analysis. 




Figure 1. dp/p for wing-body configuration and CFD A e . 


Figure 2. Ground signatures for dp/p and A e . 
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Figure 3. dp/p for wing-body configuration and A e , r . 


Figure 4. Ground signatures for dp/p and A e , r . 


The comparison shows good agreement between A e and A e?r for v < 94 (see Fig. 5). The CFD surface pressure 
distributions (see Fig. 6) do not clearly indicate why A e and A e?r become drastically different for v > 94. Note that 
the end area value of A e that is obtained from the Mach angle cuts of the CFD surface solution is completely 
dependent on the total lift, but the end area value of A e , r is affected by other three-dimensional flow properties, such 
as the pressure variations behind the aircraft. 
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IV. Matching Reversed Equivalent Area to a Target for Low-Boom Design 

The analysis in section III shows that the ground signatures for dp/p at 3BL and A e , r are nearly identical for a 
wing-body configuration (see Fig. 4). In fact, the differences between the ground signatures for dp/p at 3BL and A e?r 
are caused by two sources of numerical error — numerical error in the reverse propagation and in the conversion 
from A e , r to the near-field pressure in the forward propagation. These errors are insignificant and cause no 
significant difference between the two signatures. As a result, the problem of off-body pressure tailoring for low- 
boom design is equivalent to finding a configuration that has a low-boom A e , r - In particular, we can design a low- 
boom configuration by matching its A e?r to a low-boom target. 

The design process consists of fuselage/pod shaping for equivalent-area matching and design of experiments 
(DOE) for modification of the aft shape of A e , r - The baseline configuration was originally designed by using a 
mixed-fidelity A e matching process 2 for a cruise weight of 30,000 lb, a cruise Mach number of 1.6, and a cruise 
altitude of 45,000 ft. (See Fig. 10 in Ref. 2 for the sonic -boom analysis of this configuration.) The same mixed- 
fidelity design process that is given in Ref. 2 (except that A e is replaced with A e , r ) is used here for volume shaping to 
tailor the overall shape of A e , r , while a three-level full-factorial DOE for three twist angles of the horizontal tail is 
used to modify the aft shape of A e?r . 

To track the design progress, we used “B” to indicate the fuselage/pod shaping by BOSS 4 and “D” to indicate the 
DOE runs. With this labeling system, D3B7 is the configuration that was obtained from the third DOE run of B7, 
while B7 is the configuration that was obtained from the seventh fuselage/pod shaping iteration of the baseline 
configuration. The final design is D8B1D3B7, which has a fully shaped ground signature and a PLdB of 76.6. 

Each fuselage/pod shaping iteration that is performed with BOSS uses the following mixed-fidelity A e?r formula: 


a mixed _ ziold _i_ a volume _ a volume 
^ e,r ^e,r ' ^e^ew ^e^old 


( 12 ) 


The mixed-fidelity formula does not require the calculation of A e?r for modified configurations during the BOSS 
optimization iterations, which makes the fuselage/pod shaping process efficient. The complete inverse design 
process to obtain D8B1D3B7 can be summarized as follows: 

(a) Compute the dp/p at 3BL below the configuration. 

(b) Compute A e?r from dp/p at 3BL. 

(c) Find a low-boom equivalent-area target that is as close to A e?r as possible. 

(d) If the differences between the target and A e?r cannot be reduced by fuselage/pod volume modification, then 
go to step (g). 

(e) Modify the fuselage or the pod so that A™ xed is closer to the target configuration than A e?r (= A° ld ). 

(f) Go back to step (a). 

(g) Use DOE for three twist angles to reshape the horizontal tail. The objective function is the PLdB value of 
the ground signature for dp/p at 3BL. Pick the modified configuration that has the lowest PLdB value. 

(h) If the ground signature is not fully shaped, then go back to step (b). 
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Note that the final design D8B1D3B7 is obtained when step (h) is reached for the second time. During the 
inverse design process, the angle of attack for the Cart3D sonic-boom analysis was fixed at 0.5 deg; this value was 
determined by the baseline cruise weight of 30,000 lb. The low-boom targets were generated by a parametric 
equivalent-area target explorer. 3 




Figure 8. Comparison of ground signatures for B2 design. 



Figure 9. A e?r matching for B7 design. 




Figure 11. Comparison of A e?r for D3B7 design. 



Figure 12. Comparison of ground signatures for D3B7 design. 


After two fuselage shaping iterations, the gap between the A e?r of the configuration and the target for X < 107 was 
eliminated (see Fig. 7), which resulted in a reduction of 4 PLdB in the ground signature (see Fig. 8). Further 
fuselage shaping was performed until the aft shape of the ground signature was close to the sinusoidal shape. After 
five additional fuselage shaping iterations, the gap between the A e?r of the configuration and the target for X < 125 
was eliminated (see Fig. 9); however, the PLdB value of the ground signature for the B7 design remained higher 
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than that of the B2 design (see Fig. 10). The B7 design was accepted because the aft shape of its ground signature 
had a much lower peak overpressure and was less likely to be an A-wave than that of the B2 design. Each A e?r 
matching iteration typically required a new low-boom target. 

The differences between the A e , r of the B7 design and the target were not affected by volume shaping. Thus, a 
DOE for three twist angles of the horizontal tail was used to further modify the aft shape of the ground signature. 
The limits of the twist angles for the root, midspan, and tip airfoils were -0.2 to 0.2, -0.4 to 0.4, and -1.0 to 1.0 deg, 
respectively. A three-level full-factorial DOE was performed with 27 configuration variations. The third DOE 
configuration, with twist angles of 0.2, -0.4, and -1.0 deg for the root, midspan, and tip airfoils, respectively, was 
selected because of the aft shape of its ground signature (see Fig. 12). Figure 11 shows that A e?r was affected by the 
horizontal tail at X > 130. 




Figure 13. A e , r matching for B1D3B7 design. Figure 14. Comparison of ground signatures for B1D3B7 design. 



Figure 15. Comparison of ground signatures for D8B1D3B7 
design. 



Figure 16. Change in dp/p for D8B1D3B7 design. 


Another A e r target was generated for aft pod shaping. The pod shaping was intended to eliminate the differences 
between the A e?r of the configuration and the target for 1 17 < X < 129, but the actual CFD analysis for A e?r shows 
that the differences for 1 17 < X < 135 were only partially eliminated (see Fig. 13). The improvement in the aft shape 
of the ground signature is difficult to quantify, but a reduction of 1.5 PLdB in the perceived loudness of the ground 
signature makes the B1D3B7 design more desirable than the D3B7 design (see Fig. 14). Further fuselage and pod 
shaping did not yield any meaningful improvement of the ground signature. Another DOE for three twist angles was 
performed to tailor the aft shape of the ground signature. The same twist angle limits and DOE settings were used. 
The eighth DOE configuration, with twist angles of -0.2, 0.4, and -1.0 deg for the root, midspan, and tip airfoils, 
respectively, was selected because of the nearly perfect aft shape of its ground signature (see Fig. 15). The minor 
differences that can be seen between the dp/p distributions of the B1D3B7 and D8B1D3B7 designs demonstrate the 
sensitivity of the ground signature to the dp/p at 3BL (see Fig. 16). 

Because the angle of attack is fixed at 0.5 deg during the inverse design process, the signature for the 
D8B1D3B7 design shown in Fig. 15 is for a cruise weight of 27,800 lb. A sensitivity analysis of the ground 
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signature for the D8B1D3B7 design with respect to the cruise weight is performed for a set of angles of attack. 
Figure 17 shows that the D8B1D3B7 design still has a fully shaped ground signature and a PLdB of 78.5 at the 
desired cruise weight of 30,000 lb; Fig. 18 shows that the shapes of all A e?r are similar. These two figures show the 
sensitivity of the ground signature and the A er for the D8B1D3B7 design with respect to cruise weight changes. 



Figure 17. Comparison of ground signatures for different cruise 
weights. 



Figure 19. Surface dp/p for original and final designs. 



Figure 18. Comparison of A e?r for different cruise weights. 





Figure 21. A e?r for original and final designs. Figure 22. Ground signatures for original and final designs. 

Figures 19-22 compare the surface pressure distributions, the dp/p distributions at 3BL, the A e r values, and the 
ground signatures, respectively, for the original CFD A e design (i.e., the baseline) and the final inverse design. The 
most noticeable differences in the surface pressure distributions of the two designs are on the trailing edge of the 
wing below the nacelle and on the horizontal tail (see Fig. 19). Note that the pressure differences on the lower 
surface of the wing could be caused by the combined effect of the fuselage shaping and the difference in the angles 
of attack for the two designs. 
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V. Concluding Remarks 

The augmented Burgers equation is used to propagate the computational fluid dynamics (CFD) pressure 
distribution (dp/p) at a mid-field location to a near-field location. Then, the near-field dp/p is converted to the 
reversed equivalent area (A e , r ) by using classical far-field theory. The resulting A e , r contains much more information 
about the three-dimensional flow for sonic-boom analysis and has almost the same ground signature as the three- 
dimensional configuration when analyzed by propagation of the mid-field pressure distribution. As a result, the 
problem of tailoring a mid-field pressure distribution for a low-boom design can be solved by matching A e?r to a low- 
boom target. 

Matching the A e?r to a low-boom target has three clear advantages over tailoring a mid-field dp/p\ (1) the 
differences between the A e?r and a low-boom target indicate how far the current configuration is from a three- 
dimensional low-boom design, (2) the differences between A e?r and a low-boom target can be used to directly 
determine design regions and define design variables for low-boom shaping, and (3) finding an achievable low- 
boom target for A e?r is much easier than finding an achievable target for a mid-field dp/p. 

A low-boom supersonic demonstrator concept for a cruise Mach number of 1.6 and a cruise weight of 30,000 lb 
was developed by matching A e?r to a low-boom target. The final design has a fully shaped ground signature with a 
perceived loudness value of 78.5 PLdB. 


Appendix 

The matrices for the reversed molecular relaxation and the reversed thermoviscous absorption in the iterative 
reverse propagation process are tridiagonal and are defined as follows: 


/ l o o 

I 0 V 0v 

1-2 Sr 8r+- 


Anr(0v'£ v ) - 


Cv 2 (Ax) 


Cv 2 (Ax) 


^ 2 (At) 


\ 


1 — 2 Sr + 


2 (At ) j 


^mr(^v C,,) — 


Pc,~^ 1 ~ 2 Pc„ Pc..+ 9 


2 (At) 


2 (At) 


Pc.-- 


Cv 2 (At) 

0 


1 2 Pc, Pc, + 2 (At) 


Avr — 


/ 1 9 ' 

(p 1 — 2 (p (p 

0 


—cp 1 + 2 (p —cp 

0 ' 



> Btvr ~ 



0 

<P 1-2 <p <p 

n 't / 


0 

-<p 1 + 2 <p -<p 

n 'i 1 


0 


The constants in the above matrices are defined as 
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